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Abstract 

Exact results are obtained for random walks on finite lattice tubes with a single source and 
absorbing lattice sites at the ends. Explicit formulae are derived for the absorption probabilities at 
the ends and for the expectations that a random walk will visit a particular lattice site before being 
absorbed. Results are obtained for lattice tubes of arbitrary size and each of the regular lattice 
types; square, triangular and honeycomb. The results include an adjustable parameter to model 
the effects of strain, such as surface curvature, on the surface diffusion. Results for the triangular 
lattice tubes and the honeycomb lattice tubes model diffusion of adatoms on single walled zig-zag 
carbon nano-tubes with open ends. 
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I. INTRODUCTION 



The problem of random walks on finite lattices is fundamental to the theory of stochastic 
jrocesses jl| and has numerous applications including potential theory , electrical networks 
, atomic surface diffusion and diffusion on biological membranes • A classic problem 
in this area, which was posed by Courant et al. p in 1928, concerns random walks on finite 
planar lattices with a single source and absorbing boundaries. The exact solution for this 



problem on the square lattice was derived in 1940 jjfl. The exact solution on the triangular 
lattice was only obtained recently 0, 0] after having been considered intractable llO I. Other 
variants of the problem on the square lattice have also been solved exactly [3, 1 111 . I2I ] . 

In this paper we consider random walks from a single source on a finite lattice which 
incorporates periodic boundary conditions in one direction but absorbing boundary condi- 
tions in the other. Diffusion thus occurs on a surface with the topology of a lattice tube 
with absorbing sites at the ends. We present explicit results for each of the three regular 
lattice types; square, triangular and honeycomb. Our results include a bias parameter that 
can be adjusted away from unity to model different random walk probabilities in the cyclic 
direction around the tube compared with the axial direction along the tube. This parameter 
can be adjusted to model the effect of surface curvature, and other types of strain, on the 
surface diffusion. 

Our derivations generalize the approach developed by McCrea and Whipple [t| for planar 
square lattice random walks with absorbing boundaries. This is straightforward in the case 
of the square lattice tube however special care has to be exercised in an appropriate choice 
of co-ordinates in the case of the triangular lattice tube and the honeycomb lattice tube. 

There are two important motivations for our study. One of the motivations is to add 
to the rather small class of exactly solvable random walk lattice problems with absorbing 
boundaries, since it is still the case that: "Explicit solutions are known in only a few cases" 
1^ | . The second motivation is that the mathematics of discrete lattice diffusion problems 
may find applications in the recently realized laboratory assembly of lattice nanostructures 
(see for example, For example our problem of random walks on the triangular lattice 

tube and the honeycomb lattice tube represent models for adatom diffusion on single walled 
zig-zag carbon nano-tubes 3] with open ends - the random walks on the triangular lattice 
tube model adatom diffusion across carbon-carbon bonds and the random walks on the 
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honeycomb lattice tube model diffusion along the carbon-carbon bonds. The diffusion of 
carbon adatoms along the carbon-carbon bonds of a carbon nanotube plays a vital role in 
stabilizing and maintaining the open edge growth of nanotubes jjjjl ~v\ . Our exact results 
complement related results for diffusion on carbon nano-tubes based on i) enumeration of 

n 

random walks up to a set length [18[ and ii) microcanonical molecular dynamics simulations 

The remainder of the paper is divided into separate sections for square lattice tubes, 
triangular lattice tubes, honeycomb lattice tubes, and a section containing an example and 
discussion. 



II. SQUARE LATTICE TUBES 

Consider the standard square lattice co-ordinates (p, q) representing the intersections of 
equidistant vertical straight lines and equidistant horizontal straight lines. The expectation 
that a random walk starting from a site (a, b) visits a site (p, q), distinct from (a, b), before 
being absorbed at a finite boundary site is given by the homogeneous partial difference 
equation 

F(p,q) = j-^-[F(p + l,q)+F(p-l,q) 

+rjF(p,q+l)+riF(p,q-i)]. (1) 

The parameter 77 allows for different probabilities for walks around the tube compared with 
walks along the tube. To accommodate the source term at (a, b) we construct separate 
solutions; Fj(p,q) for q < b and Fu(p,q) for q > b. The expectation that a random 
walk starting from (a, b) visits a site {p,b), not necessarily distinct from (a,b), before being 
absorbed at a finite boundary is then given by the inhomogeneous partial difference equation 

Fj(p, b) = 5 p , a + [Fj(p +1,6) + Fj(p - 1, b) 

+ V F I (p,b-l)+ V F II (p,b + l)}. (2) 

The above difference equations are to be solved with periodic boundary conditions in the p 
co-ordinates, 

F{p,q) = F(p + m + l,q), (3) 
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absorbing boundary conditions in the q co-ordinates, 

F I (p,0) = J (4) 
F H (p,n + 1) = 0, (5) 

and matching conditions at q — b, ie., 

F I (p,b) = F II (p,b). (6) 

The method of solving inhomogeneous linear partial difference boundary value problems 
as above consists of two parts [?], 0, Q Q| • First obtain the general separation of variables 
solution to the homogeneous problem, then find an appropriate linear combination of such 
solutions to satisfy the boundary conditions and the inhomogeneous problem. A major 
difficulty in these problems can be the identification of a lattice co-ordinate system with a 
separation of variables solution that can be matched with the boundary conditions jjj. 

The homogeneous field equations for the square lattice, Eq.(JTJ), admit the separable so- 
lution 

F(p,q) = P(p)Q(q) (7) 

where 

P(p + 1) + (A - (2 + 2r]))P(p) + P(p - 1) = 0, (8) 

Q(q + l)-±Q(q) + Q(q-l) = 0, (9) 
V 

and A is the separation constant. These separated equations have general solutions 

P(p) = Aff + Bfi~ p , (10) 
Q(q) = Cu q + Du~\ (11) 

where 

<t= 2 + a,-A + yp + a? -A)»-4 | 



2r] 2 K ' 

If A = 2i] then the solutions are no longer provided by Eqs. (fTUjl and (JTTJ). In this case we 
have the solutions 

P(p) = A + Bp (14) 

Q(q) = C + Dq. (15) 



The special solutions in Eqs. (jl4|) and ()15j) do not appear in the planar lattice problems 
because they cannot satisfy absorbing boundary conditions in both p and q co-ordinates. 

With suitable linear combinations of the above solutions, Eqs. (fTU|) .(|TT )) .(fT H) .(fTo ]) . we find 
that the general solutions to the homogeneous problem that also satisfy the boundary con- 
ditions, Eqs.(J2J),(@J), and the matching condition, Eq.fjHJ), can be written in the form 

m 

Fj(p,q) = C q(b-n-l) + Y,(ake iakP + b k e- ia * p ) 

k=l 

x sinh(/3 fe ?) sinh[/3 fc (& - n - 1)] (16) 

m 

F n (p,q) = cb(q-n-l)) + Y,^ke iakP + b k e- iakP ) 

k=l 

x sinh[/3 fc (g - n - 1)] sinh(/3 fc 6) (17) 



where 



and 



2 nk 

a k = —— 18 
m + 1 



2 + 2-q = 2-q cosh f3 k + 2 cos a fc (19) 

The constants c, and are now determined by the requirement that the solutions 
satisfy the inhomogeneous equation, Eq.((2J). This step is facilitated using the identity in 
Eq. (jl9|) together with the identity in Eq. (jA.l|) and the Kronecker delta identity 



in 



£ e 27ri ^-°) fc /( m+1 ) = (m + l)5 Pja (20) 



re-arranged as 



k=0 



! I / 



,fc=l 



m + 1 2(m + 

+ ^ e -*«*P e <«*« j . (21) 

k=l J 

We thereby obtain the following solutions for the expectations that a random walk will visit 
a site (p, q) before being absorbed at an end site: 

(2 + 2r])q(n + 1-6) 



Fi(p,q) 



rj(n + l)(m + 1 
2 + 277 



— ^cos[a fc (p- a)] 

(22) 



sinh[j3 k (n + 1 — 6)] sinh(/3 fc g) 



sinh(/3 fc ) sinh[/3fe(ra + 1)] 
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FnM= (2 + a,)K» + i- 9 ) 



rj(n + l)(m + 1) 
2 + 2r] 



77(777. 



x 



— ^cos[a fc (p- a)] 
' fc=i 

sinh[/? fc (n + 1 - q)] smh((3 k b) 
sinh(/3 fc ) sinh[/3 fc (n + 1)] 

with f3 k dependent on k through Eqs.(fl8 p .([19| ) . 

The absorption probabilities G(p, q) are readily evaluated from 



It is a simple matter to show that 

n + 1 — b 



p=0 



n + 1 



and 



p=0 

so that 



m / h \ 

yj G(P , n + D = (—) 



£[G(p,0)+G(p,n + l)] = l. 

p=0 



III. TRIANGULAR LATTICE TUBES 



(23) 



G(P,0) = (24) 



To obtain the solution for triangular lattice tubes we consider the lattice co-ordinates 
(p,q) shown in Figure 1. This co-ordinate system describes two independent triangular 
lattice systems, only one of which can be accessed by a random walk from a single point 
source. In Figure 1 the sites that are accessible from the source site at (a, b) are indicated 
by filled circles. The nearest neighbour sites to the source are highlighted by open circles 
in this figure. The lattice sites that are not accessible from the source are referred to as the 
ZeI o m esh 0. The Sam e co-ordinate SySt e m n as beeo Us ed to fled ae approbate solution 
to the planar triangular lattice problem with absorbing boundaries |l0| however a different 
zig-zag co-ordinate system was required to find the exact solution |9j. In the case of the 
planar triangular lattice problem the co-ordinate system in Figure 1 allows the leakage of 
random walks from sites at p — 1 and p = m to sites at p — —1 and p = m + 2 respectively, 
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FIG. 1: Triangular lattice with (p, q) co-ordinates as used in section III. 



which are outside the absorbing boundaries at p = and p — m + 1. In the triangular lattice 
tube this leakage is prevented because the p co-ordinate is cyclic. Note that m+1 must be 
an even integer to permit periodic boundaries in the p direction. 

With the co-ordinate system shown in Figure 1 the problem on the triangular lattice tube 
is described by the homogeneous equation, for (p, q) ^ (a, b), 



F(p, q) = [F(p + 2, q) + F(p - 2, q) 



+r]F(p + 1, q + 1) + r]F(p + 1, q - 1) 
+ V F(p - 1, q + 1) + V F(p -l,g-l)], (26) 



the inhomogeneous equation, for q — b, 



Ffa, b) = 8 p , a + [F/(p + 2, b) + F T (p - 2, b) 

+r)F n (p + 1, b + 1) + T]Fi(p +1,6-1) 
+ V F H (p -1,6+1)+ r]Fi{p -1,6-1)], (27) 



and the boundary conditions, Eqs.fjHJ),®,©- 
The homogeneous equation separates as 



with solutions 



where 



P{p + 2) + \P(p + 1) - (2 + Arj)P{p) 

+XP(p-l) + P(p-2) = 0, (28) 

Q(q+l)-±Q(q) + Q(q-l) = 0, (29) 
V 



P{p) = Ae iap + B~ iap + Ce il3p + De~ iPp , (30) 
Q(q) = Ee^ + Fe™, (31) 

a = cos" 1 (-- + ^A 2 + 16 + 16r/ j , (32) 



-i (-^ + ^^ + 16 + 16^ 
/? = cos" 1 (~ - 1^2 + 16 + 16^ 



(33) 



7 = cosh-^A). (34) 

The general solution to the homogeneous problem on the triangular lattice that satisfies all 
the boundary conditions, Eqs. (J3J) , (J2J) , © , can be written as 

PiiPiO) — cq(b — n — 1) 

+dq(b — n — 1) cos[7r(p — a)] cos[7r(g — b)} 

m 1 

+J2 ( a ke iakP + b k e~ iakP ) 

k=l 

x sinh(7 fc g) sinh[7 fc (6 - n - 1)], (35) 
Fn(p,q) = cb{q - n - 1) 

+db(q — n — 1) cos[7r(p — a)] cos[7r(g — b)] 

m 1 

+Y, ^e iakP + b k e~ iakP ) 

k=l 

x sinh[7 fc (g - n - 1)] sinh(7 fc 6), (36) 

where the prime on the sum has been used to indicate that the sum does not include the 
value k = {m + l)/2 and 

"fc = rr, 37 

m + 1 
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with 

2rj cosh 7fc cos a k = 1 + 2r] — cos 2a k . (38) 
The homogeneous solutions of the form 

(A + Bq) cos[vr(p - a)] cos[n(q - b)} (39) 

in Eqs. p5jl . p5j) . are important in two fundamental ways. First they replace the null solutions 
at k = (m + l)/2 in the representation 

( afce fa*p + b k e~ iakP ) sinh( 7fc g) sinh[ 7fc (6 - n - 1)] 

and secondly they allow the appropriate zero mesh solution for lattice co-ordinates that 
cannot be accessed by a source at (a, b). 

The constants c, d, a^, are found by substituting the homogeneous solutions, 
Eqs. p5Jl .p5 j) . into the inhomogeneous equation, Eq. (j2*7j) . Using the identity in Eq.pSJ) 
together with the idendity in Eq. fjA.lj) we first obtain the intermediate result 

- (2 + Ar])5 P:a = 2rjc{n + 1) + 2r]d(n + 1) cos[vr(p - a)} 

m 1 

+2r/Y^ (a k e iakP + b k e- iakP ) cos a k sinh 7fc sinh[ 7fc (n + 1)] . (40) 

k=l 

The final result is then found by expanding the Kronecker delta as 

5 p>a = — -— + — l — cos[vr(p - a)] 
m + 1 m + 1 

1 

4 



\ [ \ ia k p -ia k a , -ia k p ia k a \ /a-, \ 



Thus we obtain 



2{m + 



F/(p , 5)= (1 + *!)«<» + 1-») 



?7(n + l)(m + 1) 
x(l + cos[7r(p — a)] cos[7r(g — 6)]) 

1 + 2^ ^ 

+— — -> cos a fc (p- a) 

7?(m+ 1)^ 

^ sinh[ 7fc (n + 1-6)] sinh( 7fc g) 
cos «fc sinh 7fc sinh[ 7fc (ra + 1)] ' 
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Fn(p,q) 



(1 + 2rj)b{n + 1 - q) 

r](n + l)(m + 1) 
x (1 + cos[7r(p — a)] cos[7r(g — b)]) 



l + 2r] 



— ^2 cos[a k (p-a)] 



k=l 



r)(m + 

sinh[7fc(ra + 1 — q)] sinh(7fc&) 



cosafc sinh7 fc sinh[7 fc (n + 1)] ' ^ ^ 

where a k and 7^ depend on k through Eqs. (p?Tj) . (|H%|) . Our results in Eqs. (J32J) , hold for 
arbitrary m except m + 1 = 0(mod 4) where singularities occur for k = ^7^, 3(»™+i)_ rj^g 
absorption probabilities G(p, q) are given by 

V 



4 ' 



G(p,0) 
G{p,n + 1) = 



2 + 4r] 



(Fj(p + 1,1) +F T (p- 1,1)), 



V 



{F n (p + l,n)+F n (p-l,n)) 



2+4:7] 

The absorption probabilities at the ends of the tube are thus 

'n + l-b s 



p=0 



n + l 



and 



with 



m / h \ 

£ G fe„ + l)=( ) 

p=0 v 7 



53[G(p,0) + G(p,n + l)] 

p=0 



IV. HONEYCOMB LATTICE TUBES 



(44) 
(45) 



Here we consider a co-ordinate system in which we label the vertices of the honeycomb 
lattice by the intersection points (p, q) of horizontal straight lines p — 0, 1, 2, . . . , m + 1 and 
vertical zig-zag lines q = 0, 1, 2, . . . , n + 1 (see Figure 2). The expectation that a random 
walk visits an interior site (p, q), distinct from the starting site (a, b), is given by the coupled 
homogeneous difference equations: 



F(p, q) 



F(p, q) 



2 + r] 



TjF(p, q + l) + F(p+ 1, q) 
+#(p-l,g)] , 

[r]F(p, q-l) + F(p+ 1, q) 
+F(p-l,q)}. 



(46) 



(47) 
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FIG. 2: Honeycomb lattice with (p, q) co-ordinates as used in section IV. 



Here F(p, q) is the expectation at sites (p, q) with nearest neighbours on the right at (p, q + 1) 
and F(p, q) is the expectation at sites (p, q) with nearest neighbours on the left at (p,q — l). 
We will refer to these distinct symmetry sites as h sites and H sites respectively. The appeal 
of this particular choice of (p, q) co-ordinates is that the difference equations for the distinct 
symmetry sites can be decoupled into separable equations for each. Indeed both F(p, q) and 
F(p, q) satisfy the same homogeneous triangular lattice equation as Eq.(J2HJ) except near the 
absorbing boundaries. This relationship between random walks on the honeycomb lattice 
and random walks on the triangular lattice is illustrated in Figure 3. In this figure the site 
at (p, q) is shown as a h site. The nearest neighbour triangular lattice sites to the site at 
(p, q) are of the same h symmetry type and co-ordinate labels for the honeycomb lattice 
and the triangular lattice match at all h sites. A similar matching occurs at H sites (as can 
be seen by inverting Figure 3). The nearest neighbours on the triangular lattice are next 
nearest neighbours on the honeycomb lattice and the probability of a random walk from 
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q-2 q-1 



q q+1 q+2 



p+2 



P+l 



♦ • f 



p 



p-i f 



P-2 



q-2 q-1 q q+1 q+2 

FIG. 3: Honeycomb lattice showing a site at (p, q) (filled box) surrounded by nearest neighbours 
on a honeycomb lattice (large circles) and nearest neighbours on a triangular lattice (small filled 
circles). The q labels at the bottom of the figure refer to the honeycomb lattice and those at the 
top refer to the triangular lattice. 



(p, q) to one of these sites is the same for random walks on the triangular lattice and the 
honeycomb lattice. 

The triangular lattice equation fails for F(p, n) due to the absorbing boundary condition 
F(p,n + 1) = and similarly the triangular lattice equation fails for F(p, 1) due to the 
absorbing boundary condition F(p, 0) = 0. To circumvent these boundary problems we use 
the homogeneous triangular lattice equation solutions for F(p, q) in the region q < b (region 
I) and we use the homogeneous triangular lattice equation solutions for F(p, q) in the region 
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q > b + 1 (region II). Thus we have homogeneous solutions of the form 



Fi{p,q) — cq + dqcos[n(p — a)]cos[7i(q — b)] 

m I 

+Y< ( a ^ iakP + he~ iakP ) sinh( 7fc g), 
fe=i 

Fn{p,q) = c(q-n-l) 

+d(q — n — 1) cos[7r(p — a)] cos[7r(g — b)] 



(48) 



k=i 



a k e 



ta k p 



b k e~ iakV ) sinh[ 7fc (g-n-l)]. 



(49) 



These homogeneous solutions satisfy the periodic boundary conditions as well as the ab- 
sorbing boundary conditions Fjj(p, n + 1) = and Fj(p, 0) = 0. 

We now consider a single point source at a type h symmetry site. Thus we have the 
inhomogeneous problem 

1 



Fj(p,b) = 5 p , a + 



T)F H (p,b+ V 



F n (p,b+1) 



2 + T] 

+F I (p + l,b)+F I (p-l,b) 

1 [r]Fj(p,b)+F n (j>+l,b + l) 



2 + T] 

+F n (p-l,b + l)]. 



(50) 



(51) 



At this stage we do not have general solutions for Fj(p, q) and F n (p, q) and so we use 
Eqs. (l46|) . (}4Tj) to write the inhomogeneous problem in the form 

Fj(p,6) = [(2 + v ) 2 5 p , a + {2 + rj)7iFn(p,b + l) 
+ jjF^p +1,6-1) + F/(p + 2, b) 
+ rjFrip - 1, b - 1) + Fjip - 2, b)] 

I [(2 + ^) 2 -2] (52) 



F H (p,b + l) = [(2 + v)vFi{p, b) 

+ 7]F n {p + 1,6 + 2) + F n {p + 2, b + 1) 

+ v F n (p - 1, b + 2) + F n (p - 2, b + 1) 

/ [(2 + ^-2]. 



(53) 
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The unknown constants c,d,a k ,b k ,c,d,a k ,b k can now be obtained by substituting the ho- 
mogeneous solutions, Eqs.()48| ) .p5 |) and the Kronecker delta identity, Eq.([41|). into the inho- 
mogeneous equations Eqs. (JH2J) . (JH3J) and equating linearly independent functions of p. The 
algebraic manipulations are simplified using the identities in Eqs. (IA.2|) . ()A.3|) . 
The resulting expressions for the expectation values are 



F/(p, q) = £±#- ( 1 - - fo + 2 ) 6 ^ q <i + caafafo - a)] cos[n(q - 6)1) 

+ 2 V {m+l) ^ {cos[a k (p~a)}smh( lk q) 

x [(?7 + 4cosh(7 fc ) cos(a fc )) sinh[7 fc (6 - n)] - 2 cos a k sinh[7 fc (6 + 1 — n)]]} 

/ {(rj cosh 7 fc cos a k - 1 - 77) (cosh [7^(72 - 2)] - cosh(7 fc n)) 

— (r] + 4 cosh7fc cosa/c) cosafc sinh7fc sinh(7fcra)} (54) 

(2 + ^)3 v^' 

+ 7TT/J {cos[a fc (p-a)]sinh[7 fc (g-n- l)]sinh(7 fc 6)} 

2r7(m+l)^ 

/ {(r/ cosh 7 fc cos a fc - 1 - r/)( cosh [7^(72 - 2)] - cosh^n)) 

— (r/ + 4 cosh7fc cosafc) cosafc sinh7fc sinh(7fcTi)} (55) 

The expectation values at type H sites in region I and h sites in region II can now be obtained 
by substituting the solutions from Eq. (j54|) into Eq. (|47j) and the solutions from Eq.(j3SJ) into 
Eq.()46|) respectively. The results are 

= aSSi) t 1 -77^72) ((-! + %- 9 )(l-«-br(p-)] cos W ,-m) 



+ 



:2 + r/) 



277(777, + 



— ^ {cos[a fc (p - a)] [77 sinh[7 fe (g - 1)] + 2 cos(a fc ) sinh(7 fc g)] 



k=l 

x [(77 + 4cosh(7 fc ) cos(a fc )) sinh[7 fc (6 - 77)] - 2 cosa fc sinh[7 fc (6 + 1 -77)]]} 
/ {(7?cosh7 fc cosa fc - 1 - 77) (cosh [7^(77 - 2)] - cosh^n)) 

— (77 + 4 cosh7fc cosafc) cosafc sinh7fc sinh(7fc77)} (56) 
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_, , N (2 + rj) 2 b [(rj + 2)(n-q) + 2 r . N1 r . ,. n 

F n (p,q) = - — L \ ^- — ' . J (l + cos 7r(p - a) cos 7r(g - 6) ) 

(2 + ^)2 ™' 
+ 2r/(m + 1) ^ { cos Mp " a )] sinh (7^) 

x (?7sinh[7 fc (g -n)] + 2 cos(a fe ) sinh[7 fc (g - n - 1)])} 

/ {(77 cosh 7^ cos a fe - 1 - 77)(cosh[7 fc (n - 2)] - cosh(7 fc n)) 

— (77 + 4 cosh7 fc cosafc) cosafc sinh7 fc sinh(7 fc n)} (57) 

The absorption probabilities are defined by 

G(p,0) = 0, (58) 

G(p,0) = -5-F 7 (p,l), (59) 
2 + 77 

G(p,n + 1) = -^-F„(p,n), (60) 
2 + 7/ 

G(p,ra+1) = 0. (61) 
The total absorption probabilities at the ends of the tubes are thus 

fa + 2)n + 2 

and 

VG(p,n + l) = fo + 2) * 
^ V ; (77 + 2)71 + 2 

p=0 

Note that the absorption probabilities in this case are functions of the bias parameter 77. The 
square lattice tube and the triangular lattice tube are symmetric with respect to left/right 
walks along the axial direction and thus the probabilities for absorption at the ends depend 
only on the initial distance from the ends at which particles are released. The honeycomb 
lattice is not symmetric with respect to left/right walks along the axial direction. As a 
consequence the absorption probabilities at the ends of the tube depend on both the initial 
distance from the ends (which also determines the symmetry type, h or H, of the initial 
lattice site) and the axial bias parameter. 



V. EXAMPLE AND DISCUSSION 



In this paper we have derived exact formulae for the expectations that a random walk 
starting at a lattice point (a, b) will visit a lattice site (p, q) on a lattice tube with absorbing 
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lattice sites on the ends. The formulae for square lattice tubes, Eqs. (j22j) . (j23j) . triangular 
lattice tubes, Eqs. p2jl . pHj) . and honeycomb lattice cubes, Eqs. (jo^j) . (l531) . (jo*DT) . (joTj) allow us 
to readily compute the expectation values for tubes of any specified size and arbitrary 
starting points. Moreover each of these solutions contains an adjustable parameter 77 that 
can be adjusted away from unity to model different random walk probabilities in the cyclic 
direction around the tube compared with the axial direction along the tube. 

As an example we consider the case of a honeycomb lattice tube with m = 17, n = 29 
and three values of 77; i) 77 = 1, ii) 77 = 1/100, and iii) 77 = 100. We have taken the source 
to be centrally located at a = 9, b = 15 in each case. The expectation values at each of the 
lattice co-ordinates have been plotted in Figure 4. As might be anticipated, for small values 
of the axial bias parameter, 77 <C 1, diffusion along the tube axis is very slow; the random 
walk cycles around the tube many times (several hundred times for 77 = 1/100) before finally 
being absorbed at one of the open ends. For large values of the axial bias parameter, 77 ^> 1, 
we might first anticipate rapid diffusion along the tube axis however the tower like plot in 
Figure 4(c) reveals that this is not the case. The random walk becomes trapped locally near 
the source as it moves back and forward between the source site and the nearest neighbour 
to the source. This effect is indeed a simple consequence of the honeycomb lattice geometry 
which has only one nearest neighbour along the axis direction. 

A further interesting calculation is the steady state profile for expectation values along 
the lattice tube after summation over p. The profile is piecewise linear with a linear increase 
from q = 1 up to q = b followed by a linear decrease from q = b to q = n. The slope of 
the linear portions is dependent on the parameter 77. Explicit expressions for this slope as 
a function of 77 can be readily evaluated from the formulae for the expectation values given 
in Eqs (54)- (57). For example for q < b we have 

m 

Ej(q) = ^2F I (p,q) + F I (p,q) 

p=0 

(V + 2) 2 / _ (n + 2)b \ _ ( V + 2) f _ ( V + 2)b \ 
77 V (v + 2)n + 2j q 2 V (v + 2)n + 2J' 

It is clear from this equation that the slope diverges as 77 — > and as 77 — > 00. For the 
example considered here with n = 29 and b = 15 we have 

(7r} + 15) (77 + 2) (2 (77 + 2)5-77) 



Ei{q) 



r/(60 + 2977) 
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FIG. 4: Expectation values for a random walk on a honeycomb tube for three different values of 
the axial bias; a) 77 = 1, b) 77 = 0.01, c) 77 = 100. 

and the slope is a minimum at 77 pa 2.035. 

The geometry of the honeycomb lattice tube that we have considered in this paper is 
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equivalent to that of a single walled zig-zag carbon nanotube with open ends. The above 
example corresponds to the (9,0) nanotube in the standard notation The bias parameter 
could thus be tuned to model the effects of strain, such as surface curvature, on diffusion 
of adatoms along the carbon-carbon bonds on zig-zag carbon nano-tubes. An interesting 
result in this connection (as shown in Figure 4(c) above) is that a random walk could be 
localized for a period of time by applying a uniform strain which favours diffusion in the 
direction of all bonds aligned with the tube axis. 

APPENDIX 

The following identites have proven useful for deriving the results in this paper 

sinh(7(6 — 1)) sinh(7(6 — n — 1)) + sinh(7&) sinh(7(6 — n)) 

= sinh(7) sinh(7(n + 1)) + 2 cosh(7) sinh(7&) sinh(7(6 — n — 1)) (A.l) 

sinh(7(6 + 1 — n)) sinh(7&) + sinh(7(6 — 1)) sinh(7(6 — n)) 

= sinh(7) sinh(7n) + 2 cosh(7) sinh(7&) sinh(7(6 — n)) (A. 2) 

sinh(7(6 — n)) sinh(7&) + sinh(7(6 — 1)) sinh(7(6 + 1 — n)) 

= (cosh(7(n - 2)) - cosh^n)) /2 (A.3) 
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